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COMPARISON  OF  NUMERICAL  AND  PHYSICAL  MODELS 
OF  WAVE  RESPONSE  IN  A  HARBOR 


PART  I :  INTRODUCTION 

Wave  Conditions  In  Harbors 


1.  Water  waves  in  harbors  are  primarily  excited  by  waves  penetrating 
through  the  harbor  entrance  or  opening.  In  engineering  practice,  waves  of 
principal  concern  for  harbor  tranquility  have  periods  in  the  range  of  about 
b  sec  to  4  min.  In  small-boat  harbors,  5-  to  20-sec  waves  are  the  primary 
concern . 

2.  The  principal  adverse  effects  of  excessive  wave  action  in  harbors 
are  associated  with  navigation  and  mooring.  Excessive  wave  action  is  caused 
by  resonance  or  insufficient  shelter  by  the  protection  breakwaters.  Resonance 
in  harbors  occurs  when  the  natural  period  of  the  harbor  is  equal  or  close  to 
an  incident  wave  period.  Resonant  oscillation  of  a  moored  vessel  occurs  when 
the  natural  period  of  the  dynamic  response  of  the  moored  vessel  is  equal  or 
close  to  an  incident  wave  period.  The  severe  wave  chop,  and  possibly  strong 
currents  associated  with  excessive  wave  action  in  harbors,  can  cause  naviga¬ 
tional  difficulty,  unacceptable  motion  of  moored  vessels,  and  undesired 
shoaling  or  erosion  within  the  harbor. 

3.  Effects  of  water  waves  in  harbors  can  be  estimated  from  field  mea¬ 
surements,  physical  model  studies,  or  mathematical/numerical  model  studies. 

All  three  sources  of  information  are  accompanied  by  assumptions  and  simplifi¬ 
cations.  Nevertheless,  it  is  accepted  that  useful  estimates  of  waves  in  an 
existing  harbor  can  be  obtained  from  field  measurements  used  in  conjunction 
with  physical  and/or  mathematical  models.  Unfortunately,  wave  measurements  in 
a  harbor  are  difficult  and  expensive  to  collect,  and  very  limited  field  data 
are  available.  For  a  proposed  harbor,  field  measurements  can  only  provide  in¬ 
cident  wave  climate;  hence,  physical  and/or  mathematical  models  must  be  used 

to  predict  wave  conditions  in  the  harbor.  The  models  are  most  useful  for  engi¬ 
neering  studies  when  they  provide  a  reasonable  representation  of  the  prototype. 

4.  Corps  of  Engineers  harbor  projects  are  designed  with  consideration 
of  the  wave  conditions  that  can  be  expected  inside  the  harbor.  Occasionally, 
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existing  harbor  projects  also  need  structural  modification  for  reasons  which 
may  include  overly  severe  wave  conditions  inside  the  harbor.  Since  costs 
associated  with  harbor  construction  and  modification  are  generally  high, 
project  engineers  need  accurate  and  efficient  models — physical  and/or 
mathematical — which  can  be  used  to  test  the  performance  of  alternative  designs 
before  construction. 


Investigation  of  Waves  In  Harbors 

5.  Historically,  both  physical  and  mathematical  models  have  been  used 
at  the  US  Army  Corps  of  Engineers  Waterways  Experiment  Station  (WES)  to  inves¬ 
tigate  waves  in  harbors.  For  example,  Bottin  (1984)  used  a  physical  model  to 
investigate  the  effects  of  proposed  harbor  improvements  on  wave  action  in  the 
small-boat  harbor  at  Barcelona,  New  York.  Houston  (1976)  used  a  numerical 
model  to  investigate  the  effects  of  proposed  improvements  on  the  response  of 
Long  Beach  Harbor  to  long-period  waves.  Many  other  studies  of  waves  in  har¬ 
bors  have  been  conducted  at  WES.  The  majority  of  these  studies  used  physical 
models  because  WES  has  excellent  laboratory  facilities  for  physical  model 
studies  and  because  numerical  models  are  relatively  new  as  practical  engineer¬ 
ing  tools.  The  advent  of  modern  computer  hardware  and  computational  methods 
has  led  to  a  dramatic  world-wide  increase  in  the  use  of  numerical  models  in 
all  fields  of  engineering.  Because  of  the  time  and  cost  efficiency  of  numeri¬ 
cal  models  for  some  studies,  their  use  is  increasing  within  the  Corps  of 
Engineers . 

6.  A  variety  of  numerical  models  has  been  used  to  investigate  waves  in 
harbors  (Lee  1969;  Chen  and  Mei  1974;  Yue,  Chen,  and  Mei  1976;  Berkhoff  1976; 
Houston  1981;  Lepelletier  1981;  Ganaba,  Welford,  and  Lee  1982;  Behrendt  and 
Jonsson  1984;  Skovgaard ,  Behrendt,  and  Jonsson  1984;  Yoshida,  Ijima,  and 
Okuzono  1984;  Matsoukis  1985;  and  Chen  1986).  These  authors  used  various 
numerical  techniques,  including  the  ray  method,  Green's  function  method, 
eigenfunction  method,  finite  difference  method,  finite  element  method,  and 
hybrid  element  method.  Recently,  WES's  Coastal  Engineering  Research  Center 
(CERC)  developed  and  implemented  the  HARBS  and  HARBD  computer  codes  to  model 
waves  in  harbors.  These  numerical  models  use  a  hybrid  element  method  to 
determine  the  wave  response  of  a  harbor  to  a  simple  time-harmonic  incident 
wave  train.  HARBS  is  used  for  shallow-water  harbors;  HARBD  is  used  for 


arbitrary-depth  harbors.  The  models  have  been  tested  against  some  laboratory 
data  (Chen  1984  and  1986)  and  have  been  used  to  predict  wave  response  in 
several  prototype  harbors  (Farrar  and  Chen  1987  and  Mathiesen  1987). 

7.  The  Shore  Protection  Manual  (SPM)  (1984)  presents  diffraction  dia¬ 
grams  and  methods  of  implementing  them  which  can  be  used  to  study  wave  propa¬ 
gation  into  a  harbor  if  the  problem  description  meets  a  variety  of  restrictive 
criteria.  The  methods  are  useful  only  for  the  simplest  problems.  It  was 
realized  that  diffraction  diagrams  were  not  strictly  applicable  to  the  problem 
of  determining  wave  conditions  in  Barcelona  Harbor.  Nevertheless,  it  was 
deemed  worthwhile  to  determine  wave  conditions  inside  the  harbor  by  using  dif¬ 
fraction  diagrams  and  then  to  compare  the  results  with  results  from  the  physi¬ 
cal  and  numerical  models.  The  method,  results,  and  comparisons  are  presented 
in  Appendix  A. 


8.  The  purpose  of  the  study  herein  is  to  determine  the  wave  response  at 
Barcelona  Harbor,  New  York,  using  the  HARBD  numerical  model  and  then  to  com¬ 
pare  the  numerical  model  results  with  data  from  the  physical  model  study 
described  by  Bottin  (1984).  This  study  is  intended  to  provide  technical 
information  which  will  help  establish  preliminary  engineering  guidelines  for 
conducting  harbor  wave  studies;  it  is  not  intended  to  provide  information  con¬ 
cerning  the  time  and  cost  requirements  of  the  two  modeling  approaches. 
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PART  II:  PHYSICAL  ANT)  NUMERICAL  MODEL  STUDIES 
OF  BARCELONA  HARBOR,  NEW  YORK 


Physical  Model  Study 

9.  At  the  request  of  US  Army  Engineer  District,  Buffalo  (NCB) ,  a 
physical  model  study  of  wave  conditions  in  Barcelona  Harbor,  New  York,  was 
conducted  at  WES  during  late  1983  and  early  1984.  The  effects  of  a  variety  of 
proposed  harbor  improvements  on  wave  action  in  the  harbor  were  investigated. 

10.  The  model  study  is  described  in  detail  by  Bottin  (1984).  A  brief 
description  of  relevant  aspects  of  the  study  is  given  here.  The  model  was  a 
I:60-scale  undistorted  model.  A  time  scale  of  1:7.75  was  required  by  the 
Froude  model  law.  Figure  1  shows  the  model  layout.  Model  tests  were  per¬ 
formed  for  60  different  harbor  configurations,  three  of  which  were  considered 
as  part  of  the  numerical  model  study  described  in  paragraphs  28  through  38. 
Only  waves  with  periods  between  4.9  and  10.1  sec  were  considered. 

11.  Model  tests  were  conducted  for  monochromatic  and  irregular  waves. 
During  development  and  selection  of  an  optimum  plan,  monochromatic  model  tests 
were  made  for  each  harbor  configuration  using  a  variety  of  water  levels  and 
incident  waves  with  different  heights,  directions,  and  periods.  Irregular 
wave  tests  were  made  only  for  the  optimal  proposed  plan  (Plan  58  in  Bottin' s 
(1984)  report).  The  irregular  wave  tests  represented  CERC's  initial  develop¬ 
ment  of  the  capability  to  generate  irregular  wave  conditions  in  a  three- 
dimensional  model.  (The  capability  to  conduct  irregular  wave  tests  in  wave 
flumes  had  previously  been  developed.)  For  each  test,  wave  records  were  taken 
at  14  gages  within  or  near  the  harbor.  The  locations  of  Gages  3  through  14 
were  the  same  for  all  harbor  configurations.  The  locations  of  Gages  1  and  2 
varied  slightly  for  different  harbor  configurations,  but  they  were  always  near 
the  harbor  entrance.  The  locations  of  the  gages  are  shown  in  Figures  2 
through  4. 

12.  For  the  monochromatic  tests,  most  wave  records  had  some  variation 
of  individual  wave  heights  and  periods.  For  each  test,  a  crest-to-t rough  wave 
height  analysis  technique  (Turner  and  Durham  1984)  was  used  to  calculate 

,  the  average  of  the  one-third  highest  waves  in  the  record,  at  each  gage. 
The  technique  closely  duplicates  analysis  techniques  used  for  monochromatic 
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attenuation  caused  by  viscous  bottom  friction  in  the  model.  The  adjustment, 
which  is  a  function  of  temperature,  is  estimated  to  be  in  the  range  of  4  to 


19  percent.  Values  of  the  adjusted  are  used  to  characterize  the  wave 


conditions  at  each  gage  (Bottin  1984).  The  use  of  the  adjusted  H^^  to 


represent  the  wave  height  at  a  gage  in  the  physical  model  study  is  regarded  as 
a  good  estimate  for  engineering  purposes  such  as  plan  evaluation.  However, 
the  adjusted  H^^  does  not  fully  characterize  the  model  data.  (See 
paragraphs  19  through  27.) 

13.  Scale  effects  (Hudson  et  al.  1979)  in  a  small-scale  physical  model 
are  not  limited  to  viscous  energy  dissipation  at  the  bottom.  Other  scale  ef¬ 
fects  are  related  to  transmission  through  and  reflection  from  breakwaters,  re¬ 
flection  from  the  sloped  transition  section  between  the  generator  pit  and  the 
modeled  bathymetry,  and  reflections  from  the  wave  generator  and  model  basin 
boundaries.  Physical  models  are  carefully  designed  to  minimize  these  scale  ef¬ 
fects.  Operationally,  monochromatic  wave  test  data  are  acquired  and  analyzed 
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prior  to  the  time  at  which  waves  reflected  from  the  model  boundary  or  wave 
generator  return  to  the  harbor  area. 

14.  Wave  generator  characteristics  also  affect  monochromatic  and  irreg¬ 

ular  wave  generation.  In  practice,  as  the  stroke  of  the  wave  generator  in¬ 
creases,  wave  generators  designed  to  produce  monochromatic  waves  usually  tend 
to  generate  energy  at  harmonics  of  the  fundamental  frequency.  The  wave  record 
from  a  gage  located  near  the  center  and  10  ft*  in  front  of  the  wave  generator 
for  a  6.7-sec  (prototype)  wave  is  shown  in  Figure  5a.  A  0.2-in.  (model) 
stroke  was  used  to  generate  this  wave.  The  H ^ ^  of  this  wave  is  1.66  ft 
(prototype!,  and  the  mean  amplitude  is  1.64  ft  (prototype).  The  corresponding 
energy  spectrum,  shown  in  Figure  5b,  indicates  very  little  wave  energy  at 
harmonic  frequencies.  A  taper  function  was  applied  to  the  data  during  the 
frequency  analysis,  but  smoothing  techniques  were  not  used.  The  wave  record 
from  the  same  gage  for  a  larger  amplitude  wave  having  a  period  of  5.3  sec 
(prototype)  is  shown  in  Figure  6a.  A  0.5-in.  (model)  stroke  was  used  to  gen¬ 
erate  this  wave.  The  for  this  wave  is  6.74  ft  (prototype),  and  the 

mean  wave  height  is  6.67  ft  (prototype).  Generation  of  a  small  amount  of 
energy  at  the  first  harmonic  is  evident  in  the  corresponding  energy  spectrum 
shown  in  Figure  6b.  Spectral  analysis  results  also  indicate  wave  energy  at 
frequencies  adjacent  to  the  fundamental  frequency.  The  amount  of  energy  at 
adjacent  frequencies  is  relatively  small.  However,  the  presence  of  even  a 
relatively  small  amount  of  wave  energy  away  from  the  intended  frequency  can 
significantly  affect  wave  propagation  and  transformation.  It  should  also  be 
noted  that  energy  transfer  to  harmonics  may  occur  in  shallow  water  in  both  the 
physical  model  and  prototype. 

15.  If  paddle  wave  generators  are  used  for  irregular  wave  tests,  incor¬ 
rect  boundary  conditions  may  cause  spurious  long  waves  (Jensen  and  Warren 
1986).  Depending  on  the  test  conditions,  the  total  long-period  wave  energy  in 
a  model  may  be  increased  or  decreased.  Long-period  wave  energy  associated 
with  wave  setdown  under  wave  groups  will  occur  and  should  be  expected  during 
irregular  wave  tests.  As  noted  by  Jensen  and  Warren,  long-period  waves  will 
be  reflected  and  re-reflected  from  nondissipating  boundaries  such  as  the  wave 
generator  and  basin  boundaries,  but  limited  comparisons  with  prototype  data 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


seem  to  indicate  long-wave  energy  in  physical  models  can  be  simulated  in  a 
realistic  way.  Additional  studies  are  needed  in  this  area.  Theoretical  re¬ 
sults  to  date  consider  piston/paddle  wave  generators  and  not  the  plunger  gen¬ 
erator  used  for  the  Barcelona  Harbor  study. 

16.  Most  of  the  wave  records  in  the  vicinity  of  the  harbor  from  the 

monochromatic  model  tests  are  not  as  regular  cs  the  records  from  deeper  water 
near  the  wave  generator.  Instead,  they  show  the  effects  of  nonlinearity  and 
exhibit  some  degree  of  wave  grouping.  Some  wave  records  also  contain  small- 
amplitude,  high-frequency,  random  oscillations  which  result  from  wave  breaking 
and  overtopping  of  the  breakwaters.  The  determined  from  the  analysis 

described  in  Paragraph  12  does  not  fully  characterize  these  wave  records  and 
was  not,  by  itself,  adequate  for  the  purpose  of  making  comparisons  between 
physical  and  numerical  model  results. 

17.  In  Bottin's  report  only  H^^  values  were  presented.  For  this 
study  the  physical  model  data  were  reanalyzed  using  a  slightly  different  proce¬ 
dure  which  permitted  investigation  of  the  full  distribution  of  wave  heights. 
Wave  heights  were  determined  by  performing  a  simple  upcrossing  analysis:  two 
consecutive  upcrossings  delineate  a  wave,  the  height  of  which  is  equal  to  the 
difference  in  the  maximum  and  minimum  water  surface  elevations  between  the 


upcrossings.  From  the  resulting  wave  height  distributions,  the  maximum  wave 

height  H  ,  the  minimum  wave  height  H  .  ,  the  wave  height  with  an  exceed- 

max  min 

ance  probability  of  50  percent  H  ,  and  the  average  of  the  one-third 

highest  waves  were  determined.  These  parameters  were  used  to  charac¬ 

terize  the  oscillations  at  each  wave  gage.  Because  of  the  limitation  of  time 
and  resources,  no  effort  was  made  to  conduct  further  statistical  analyses  of 
the  wave  data. 

18.  The  resulting  H,  ,  values  (presented  later  in  this  report)  are 

1/3 

generally  lower  than  those  presented  by  Bottin  (1984)  because  the  simple 
upcrossing  analysis  used  for  this  study  sometimes  classified  small-amplitude 
oscillations  between  crest  and  trough  as  individual  waves.  Furthermore,  the 
physical  model  wave  heights  presented  in  this  report  have  not  been  adjusted 
using  Keulegan's  equation.  For  comparison  purposes,  the  adjustment  was  not 
necessary  because  bottom  friction  coefficients  used  in  the  numerical  model 
were  chosen  which  were  appropriate  for  the  scale  of  the  physical  model.  In 


general,  from  Bottin  (1984) 

mined  for  this  study.  This  result 


lies  between  the  H  and  H,  deter- 

max  1/3 

is  a  matter  of  convention  and  is  not 


important  relative  to  the  comparisons  of  physical  and  numerical  model  results 
presented  in  Part  III. 

19.  Typical  wave  records  from  inside  the  harbor  are  shown  in 
Figures  7(a)  through  10(a).  For  comparison,  Figure  11a  is  a  typical  wave 
record  from  deeper  water  near  the  wave  generator  for  a  relatively  low  ampli¬ 
tude  wave.  The  corresponding  cumulative  probability  distributions  (CPD's)  of 

i  wave  heights,  derived  from  the  results  of  simple  upcrossing  analyses,  are 

shown  in  Figures  7b  through  lib.  The  variation  of  wave  heights  is  easily  seen 
in  the  CPD's. 

20.  For  each  run,  the  values  of  H  ,  H  .  ,  H~  ,  and  H,  ,»  , 

max  mm  0.50  1/3 

were  normalized  by  the  value  of  the  generated  wave  height  of  the  run,  as 

listed  in  Table  l.  The  results  are  plotted  along  vertical  bars  in  Figures  7c 

through  He.  Hereafter,  references  to  H  ,  H  .  ,  H„  ,  and  H,  are 

max  mm  0.50  1/3 

to  cheir  normalized  values.  Essentially,  the  vertical  bar  plots  are  simple, 
graphical  representations  of  the  CPD’s;  they  provide  a  convenient,  concise 
means  of  reporting  wave  height  data  from  the  monochromatic  physical  model 
tests.  Similar  plots  are  used  in  Part  III  to  make  comparisons  between  the 
results  of  monochromatic  physical  model  tests  and  the  numerical  model. 

21.  The  positions  of  the  four  wave  heights  on  the  vertical  bar  plots 

depend  on  the  nature  of  the  corresponding  wave  records.  For  a  wave  record 

which  is  a  simple  sinusoid,  the  CPD  is  a  step  function;  and  H  ,  H  ,  , 

max  min 

Hq  ,_q  ,  and  plot  at  exactly  the  same  point  on  the  vertical  bar  plot. 

However,  none  of  the  wave  records  considered  in  this  report  plot  as  simple 
sinusoids.  Deviations  from  the  simple  sinusoid  can  generally  be  explained  in 
terms  of  wave  grouping,  nonlinearity,  and  the  small-amplitude,  high-frequency, 
random  oscillations  caused  by  wave  breaking  and  overtopping. 

22.  Wave  records  in  and  near  the  harbor  from  the  monochromatic  physical 
model  tents  can  be  roughly  sorted  into  the  four  cases  typified  by  Figures  7 
through  10;  however,  Figures  8  and  9  typify  the  majority  of  wave  records.  In 
paragraphs  23  through  26  the  four  cases  are  discussed.  Attention  is  given  to 
the  manner  in  which  nonlinearity,  wave  grouping,  breaking,  and  overtopping 
manifest  themselves  in  the  CPD's  and  the  vertical  bar  plots.  The  discussion 
is  given  to  provide  background  information  which  will  allow  the  reader  to  more 
easily  understand  the  comparisons  between  the  physical  and  numerical  model 
results  which  are  presented  in  Part  ill. 

23.  The  wave  record  in  Figure  7  is  not  quite  a  simple  sinusoid. 
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Cumulative  Probability 


a.  Wave  record 


Table  1 

Characteristics  of  15  Monochromatic  Runs* 


Water 

Wave 

Wave 

Photo  No. 

Run 

Harbor  Wave 

Level 

Period 

Height 

in  Bottin's 

No. 

Configuration  Direction 

ft 

sec 

ft 

(1984)  Report 

I 

Base  Test  2  West 

+5.5 

9.9 

6.2 

NA 

2 

13.9 

26 

3 

North 

+4.0 

7.5 

4.5 

NA 

4 

9.6 

38 

5 

Northeast 

+3.0 

6.7 

2.9 

NA 

6 

5.7 

44 

7 

Plan  42  West 

+5.5 

9.9 

6.2 

NA 

8 

13.9 

96 

9 

North 

+4.0 

7.5 

4.5 

NA 

10 

9.6 

102 

1  1 

Northeast 

+3.0 

6.7 

2.9 

NA 

12 

5.7 

105 

13 

Plan  58  West 

+6 . 5 

7.7 

7.9 

131 

14 

North 

+4.0 

7.5 

9.6 

132 

15 

Northeast 

+3.0 

6.7 

5.7 

133 

*  Wave 

characteristics  are  values  for  the 

prototype 

waves . 

Instead 

,  it  exhibits  the  effects  of 

nonlinearity  (as 

indicated  by 

the  steep 

peaks  and  by  the  small  perturbations 

near  the  mean  water  level)  and  shows  a 

slight 

degree  of  wave  grouping.  The 

CPD  is 

very  steep,  and  the  four  wave 

heights 

are  closely  spaced  on  the  vertical 

bar  plot . 

In  this  example,  though 

not  for 

all  cases,  the  upcrossing  analysis 

has  not  identified  any 

of  the  small 

perturbations  as  a  wave.  Every  wave 

identified  by  the  upcrossing 

analysis  has 

a  period  which  is  essentially  equal 

to  the 

period  of 

the  incident 

wave .  The 

apparent  grouping  causes  the  CPD  to 

deviate 

slightly 

from  a  step 

function  and 

the  vertical  bar  not  to  collapse  to 

a  singl 

e  point. 

If  the  upcrossing  analy- 

sis  had  identified  one  of  the  small  perturbations  as  an  individual  wave,  then 

H  .  would  have  plotted  much  lower  on  the  vertical  bar  plot, 
min  r  r 

24.  The  wave  record  in  Figure  8  is  similar  to  the  one  in  Figure  7  but 
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exhibits  a  more  pronounced  beat  pattern.  As  a  result  the  CPD  is  less  steep, 
and  the  four  wave  heights  on  the  vertical  bar  plot  are  more  widely  spaced. 

25.  The  wave  record  in  Figure  9  is  similar  to  the  one  in  Figure  8  but 

is  more  sharply  peaked,  having  some  small-amplitude  high-frequency  random 

oscillations .  The  CPD  is  much  less  steep  and  shows  two  significant  segments, 

one  below  H/H  =  0.5  and  the  other  above.  (H  is  the  root  mean  square 
rms  rms 

wave  height.)  Wave  heights  below  H/H  =  0.5  primarily  represent  the 

rms 

small-amplitude  high-frequency  random  oscillations.  The  variation  of  wave 

heights  above  H/H  =  0.5  results  from  the  beat  pattern.  On  the  vertical 
rms  r 

bar  plot,  H  ,  H,  ,  and  H-  r  are  near  one  another:  and  their  separa- 
r  max  1/3  0.5  1 

tion  gives  a  measure  of  the  degree  of  development  of  the  beat  pattern.  H  .  , 

min 

which  represents  the  small-amplitude  oscillations,  is  much  lower. 

26.  In  the  wave  records  shown  in  Figures  7a  through  9a,  waves  with  fre¬ 
quencies  equal  to  the  generated  wave  frequency  are  apparent.  Also  apparent  is 
that  the  majority  of  the  wave  energy  is  associated  with  the  frequency  of  the 
generated  wave.  The  wave  record  in  Figure  10,  however,  appears  to  be  domi¬ 
nated  by  small-amplitude,  high-frequency,  random  oscillations.  The  CPD  shows 
the  large  number  of  small-amplitude  high-frequency  random  oscillations  picked 
up  by  the  upcrossing  analysis.  The  four  wave  heights  on  the  vertical  bar  plot 
are  widely  and  evenly  spaced  indicating  the  irregular  nature  of  the  wave 
record . 

27.  Energy  spectra  for  each  of  the  four  wave  records  in  Figures  7a 
through  10a  are  shown  in  Figures  12  through  15,  respectively.  Results  for  the 
first  three  wave  records  show  that  the  wave  energy  is  concentrated  at  the  fre¬ 
quency  of  the  incident  wave  with  much  less  energy  at  harmonics.  Analysis  for 
the  wave  record  in  Figure  10  shows  a  similar  trend,  though  some  energy  is 
scattered  throughout  the  spectrum.  Hence,  even  when  breaking  and  overtopping 
are  pronounced,  much  of  the  energy  remains  at  the  frequency  of  the  generated 
wave.  Spectral  analyses  show  the  transfer  of  energy  to  harmonic  frequencies 
and  the  complexity  of  wave  pattern^  which  resulted  even  when  the  generated 
wave  was  very  nearly  monochromatic. 

Numerical  Model  Study 

28.  The  computer  code  HARBD  was  used  to  model  the  wave  response  of 
Barcelona  Harbor  to  a  varietv  of  incident  wave  conditions.  HARBD  is  a  hvbrid 


Spectral  Density  (ftJ/Hz) 


frequency  (Hz) 


Figure  14.  Energy  spectrum  for  the 
wave  record  in  Figure  9 


frequency  (Hz) 

Figure  15.  Energy  spectrum  for  the  wave 
record  in  Figure  10 
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element  model  for  calculating  time-harmonic  linear  wave  oscillations  in  and 
near  coastal  harbors  of  arbitrary  water  depth.  The  model  includes  the  effects 
of  bottom  friction  and  boundary  absorption.  The  bottom  friction  is  assumed  to 
be  proportional  to  flow  velocity  with  a  phase  difference.  The  boundary  ab¬ 
sorption  is  formulated  similar  to  the  impedance  condition  in  acoustics  and  is 
expressed  in  terms  of  wave  number  and  reflection  coefficient  of  the  boundary. 
In  the  model  the  entire  water  domain  is  divided  into  near  and  far  regions  by 
an  artificial  semicircular  boundary.  The  near  region  encompasses  the  harbor 
and  the  concerned  marine  structure  and  bathymetry.  The  far  region  is  assumed 
to  have  a  straight  coastline  and  constant  water  depth.  A  variational  princi¬ 
ple  with  a  proper  functional  is  established  to  ensure  the  matching  conditions 
are  satisfied  along  the  artificial  semicircular  boundary.  A  conventional 
finite  element  approximation  is  used  for  the  solution  in  the  near  region,  and 
an  analytical  solution  is  used  in  the  far  region. 

29.  HARBD,  and  similar  models  using  the  same  numerical  technique,  have 
been  tested  against  analytical  solutions  and  laboratory  data  (Chen  and  Mei 
1974  and  Chen  1986).  The  results  are  excellent,  but  all  the  comparisons  were 
made  for  simple  bathymetry  and  within  the  regime  of  linear  wave  theory.  As  an 
engineering  tool,  the  model  has  been  used  to  predict  wave  response  in  several 
harbors  (Bottin,  Sargent,  and  Mize  1985,  Farrar  and  Chen  1987,  and  Mathiesen 
1987).  The  predictions  are  reasonable,  but  no  prototype  data  are  available 
for  comparisons. 

30.  The  HARBD  model  is  intended  to  model  harbor  waves  which  can  be  rea¬ 
sonably  presumed  to  be  governed  by  the  mild  slope  equation.  Therefore,  as  the 
waves  become  increasingly  nonlinear  and/or  gradients  in  the  bathymetry  become 
very  steep,  the  accuracy  of  the  model  is  expected  to  decrease.  Fortunately, 
many  harbor  wave  response  problems  are  included  in  the  class  of  problems  for 
which  HARBD  is  applicable. 

31.  Another  limitation  of  the  HARBD  model  is  that  the  near-region 
finite  element  grid  can  only  cover  a  finite  extent  of  the  open  water  surround¬ 
ing  the  harbor.  This  situation  is  analogous  to  that  of  the  finite  basin  size 
limitation  of  physical  models.  The  limitation  results  because  the  prototype 
bathymetry  and  shoreline  configuration  cannot  be  modelled  exactly  in  the  far 
region  (see  paragraph  27).  Furthermore,  some  preliminary  HARBD  model  runs 
indicated  that  wave  heights  inside  the  harbor  were  slightly  sensitive  to  the 
location  of  the  artificial  semicircular  boundary.  Theoretically,  the 


limitation  could  be  removed  by  placing  the  artificial  semicircular  boundary  so 
that  it  encloses  all  important  bathymetry  and  shoreline,  thus  increasing  the 
number  of  grid  cells  in  the  near  region.  The  limitation  imposes  itself  only 
when  the  capacity  of  the  available  computer  is  insufficient  and/or  the  cost  of 
data  management  and  computations  is  too  great.  Grid  size  and  spacing  can  be 
chosen  which  will  allow  accurate  and  economical  solutions  for  many  problems. 

32.  A  possible  difficulty  in  using  HARBD  is  that  bottom  friction  and 
reflection  coefficients  must  be  specified.  Presumably  reasonable  choices  for 
reflection  coefficients  can  be  made  based  on  information  in  the  SPM  (1984)  or 
based  on  experience.  The  sensitivity  of  the  model's  accuracy  to  the  choice  of 
bottom  friction  and  reflection  coefficients  depends  on  the  application.  For 
example,  solutions  for  wind  waves  in  relatively  deep  water  are  extremely  in¬ 
sensitive  to  the  bottom  friction  coefficients,  whereas  solutions  for  the 
amplifications  of  the  resonant  modes  of  a  harbor  are  at  least  slightly  sensi¬ 
tive  to  the  bottom  friction  coefficients. 

33.  To  obtain  a  reasonable  variety  of  tests  using  the  HARBD  model, 
fifteen  monochromatic  runs  and  one  irregular  wave  run  were  made.  Wave  charac¬ 
teristics  and  harbor  configurations  were  chosen  from  the  total  set  of  runs 
made  during  the  physical  model  study.  The  three  harbor  configurations  consid¬ 
ered  (Base  Test  2,  Plan  42,  and  Plan  58  in  Bottin's  (1984)  report)  are  shown 
in  Figures  2  through  4,  respectively.  Table  1  lists  the  harbor  configuration, 
water  level,  and  incident  wave  period,  height,  and  direction  of  each  of  the 

15  monochromatic  runs’.  Tabulated  water  levels  are  referred  to  the  same  lwd 
used  in  the  physical  model  study. 

34.  Finite  element  grids  representing  each  harbor  configuration  are 
shown  in  Figures  16  through  18.  The  location  of  the  artificial  semicircular 
boundary  is  the  same  for  all  three  finite  element  grids  and  is  shown  super¬ 
posed  on  the  model  layout  in  Figure  19.  Notably,  the  finite  element  grid  does 
not  cover  the  entire  modeled  bathymetry  of  the  physical  model  because  wave 
transformation — principally  refraction  and  shoaling — which  occurs  before  the 
waves  reach  the  finite  element  grid  was  not  simulated  by  the  numerical  model. 
Ideally,  the  grid  would  have  included  the  entire  bathymetry  of  the  physical 
model.  Unfortunately,  at  the  time  of  the  study,  the  cost  of  data  management 
and  computing  for  a  larger  grid  exceeded  the  resources  of  the  study.  This 
problem  has  been  discussed  for  general  applications  of  HARBD  (paragraph  30). 

To  account  for  shoaling,  shoaling  coefficients  were  calculated  for  each  wave 
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Figure  18.  Finite  element  grid  for  Plan  58  harbor  configuration 


period  based  on  the  change  in  depth  from  the  wave  generator  pit  to  the  edge  of 
the  finite  element  grid.  Shoaling  coefficients  for  the  9.9-sec,  7.5-sec,  and 
6.7-sec  waves  were  1.14,  1.07,  and  1.03,  respectively.  Wave  heights  calcu¬ 
lated  by  HARBD  based  on  the  incident  wave  heights  in  Table  1  were  adjusted  by 
multiplving  by  the  appropriate  shoaling  coefficient.  No  attempt  was  made  to 
account  for  refraction  between  the  generator  pit  and  the  edge  of  the  grid; 
hence,  wave  transformation  between  the  generator  pit  and  edge  of  the  finite 
element  grid  is  a  cause  of  discrepancies  between  the  physical  and  numerical 
model  results. 

35.  The  total  numbers  of  nodes  and  elements  for  each  harbor  configura¬ 
tion  are  6,786  and  13,108,  6,693  and  12,877,  and  6,704  and  12,888,  for  Base 
Test  2,  Plan  42,  and  Plan  58,  respectively.  The  grids  provided  a  resolution 
of  approximately  5  grid  cells  per  wavelength  for  the  shortest-period  (6. 7-sec) 
wave  considered.  The  approximate  resolution  for  the  7.5-sec  and  9.9-sec  waves 
was  6  and  8  grid  cells  per  wavelength,  respectively.  As  discussed  in  Part  IV, 
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MOOEL  LIMITS 


semicircular  boundary  of  the  HARBD  finite  element  grid 


the  results  of  this  study  suggest  a  minimum  resolution  of  7  grid  cells  per 
wavelength  to  obtain  the  best  results  from  HAFBD. 

36.  The  bathymetry  was  obtained  from  the  same  hydrographic  surveys  used 
to  design  the  physical  model.  The  elevation  of  the  lake  floor  specified  for 
the  far  region  was  -15  ft.  The  seawall  around  the  city  dock,  along  the 
improved  western  area  of  the  harbor,  and  along  the  original  breakwaters  was 
assigned  a  reflection  coefficient  of  0.95.  The  rubble-mound  breakwater  exten¬ 
sions  and  wave  absorbers  were  assigned  a  reflection  coefficient  of  0.40. 
Elsewhere,  along  the  unimproved  shore,  a  reflection  coefficient  of  0.0  was 
assigned  because  waves  are  essentially  broken  and  absorbed  there.  The  bottom 


friction  coefficient  defined  by  Chen  and  Houston  (1987)  was  assigned  to  be 
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0.1.  The  values  of  the  reflection  and  bottom  friction  coefficients  are 
typical  of  prototype  situations  and  are  thought  to  be  reasonable  choices  for 
the  present  study.  However,  the  sensitivity  of  the  solutions  to  changes 
in  these  coefficients  has  not  been  investigated. 

37.  This  study  is  the  first  adaptation  of  HARBD  to  spectral  modeling. 
One  irregular  wave  run  was  made.  The  Plan  58  harbor  configuration  was  used. 
Wave  approach  was  from  the  northwest,  and  the  still-water  level  was  at  an  el 
of  +5.0  ft.  The  energy  density  spectrum  generated  by  the  wave  generator  dur¬ 
ing  physical  model  testing  was  used  as  input  to  the  HARBD  numerical  model. 

The  spectrum  is  shown  in  Figure  20.  To  simulate  this  run  using  HARBD,  the 
spectrum  was  divided  into  5  frequency  bands  centered  at  0.116  Hz,  0.138  Hz, 
0.161  Hz,  0.187  Hz,  and  0.211  Hz,  as  shown  in  Figure  20.  Higher  frequencies 
were  not  included  because  they  are  less  significant  in  most  engineering  prac¬ 
tice  and  because  the  finite  element  grids  shown  in  Figures  16  through  18  were 
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Figure  20. 
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prepared  for  wave  frequencies  up  to  about  0.211  Hz;  they  could  not  provide  the 
requisite  resolution  for  investigating  higher  frequencies.  Indeed,  resolution 
of  the  two  highest  frequencies  is  less  than  that  suggested  by  results  of  this 
study  (see  Part  IV).  The  decision  to  use  5  frequency  bands  and  the  choice  of 
frequencies  was  completely  arbitrary.  The  rationale  was  that  the  5  selected 
frequencies  provide  a  reasonable  representation  of  the  range  of  frequencies  in 
which  significant  energy  is  present.  HARBD  was  run  for  each  of  the  five  fre¬ 
quencies,  and  the  wave  amplitudes  at  each  gage  were  determined.  At  each  gage, 
the  energy  density  corresponding  to  each  of  the  five  frequencies  was  calcu¬ 
lated  as  the  energy  density  of  the  generated  spectrum  at  that  frequency  times 
the  square  of  the  amplification  factor  for  the  gage.  The  amplification  factor 
at  a  gage  is  defined  as  the  ratio  of  the  wave  height  at  the  gage  to  the  inci¬ 
dent  wave  height. 

38.  HARBD  was  run  on  a  CDC  Cyber  Model  205  series  501  computer.  The 
central  processing  unit  (CPU)  time  for  each  monochromatic  run  was  about  2  min. 
Each  run  required  approximately  700  system  billing  units  (SBU)  which,  at  the 

time  of  the  study,  translated  to  a  cost  per  run  of  $63.70,  assuming  a  P2  (i.e. 

overnight)  job  priority.  This  information  is  given  only  to  provide  a  rough 
indication  of  the  cost  of  a  single  HARBD  model  run.  No  inferences  concerning 
the  total  cost  of  a  study  can  be  made  based  solely  on  this  information.  The 

total  cost  of  a  study  depends  on  the  nature  of  the  study  and  includes  the  cost 

of  preparing  the  finite  element  grid(s),  debugging  and  managing  the  input  and 
output  data,  and  managing  the  study  as  a  whole.  The  cost  of  grid  generation 
and  data  management  depends  on  the  efficiency  of  the  preprocessing  and  post¬ 
processing  software  which  are  areas  of  active  development  at  the  present  time. 
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PART  III:  RESULTS  AND  COMPARISONS 


39.  The  fifteen  monochromatic  wave  runs  listed  in  Table  1  and  one  ir¬ 
regular  wave  run  were  made  using  the  HARBD  numerical  model.  Wave  response, 
which  includes  the  amplification  factor  and  phase  at  each  node  on  the  finite 
element  grid,  was  computed.  In  this  Part,  the  HARBD  results  are  presented  and 
compared  with  the  results  of  the  physical  model. 

Monochromatic  Waves 


40.  In  a  qualitative  sense,  the  major  features  of  the  wave  response  in 
and  around  the  harbor  were  similarly  predicted  by  the  physical  and  numerical 
models  for  most  runs,  even  when  breaking  and  overtopping  in  the  physical  model 
were  pronounced.  A  typical  example  is  shown  in  Figure  21  which  represents  a 


snapshot  of  the  water  surface  elevation  at  time  t  =  0  for  Run  15.  This  ex¬ 
ample  corresponds  to  Photo  1  (Photo  133,  Bottin  (1984)).  In  the  plot,  the 
isopleths  are  lines  of  equal  water  surface  elevation.  The  values  given  are 
the  water  surface  elevation  above  the  still-water  level  divided  by  the  ampli¬ 
tude  of  the  generated  wave.  Both  the  numerical  and  physical  models  indicate 
that  waves  outside  the  harbor  are  generally  high  and  that  there  is  a  complex 
interference  pattern  near  the  breakwaters.  At  the  harbor  entrance,  wave 
heights  are  about  one-third  as  high  as  outside  the  harbor.  Wave  heights  in¬ 
side  the  harbor  are  significantly  less  than  outside,  and  the  harbor  is  rela¬ 
tively  calm. 

41.  Strictly  speaking,  straightforward  quantitative  comparisons  could 

not  be  made  because  of  the  irregularity  of  the  physical  model  wave  records. 

Still  it  was  possible  to  make  meaningful  comparisons  at  each  gage  for  each 

run.  The  comparisons  are  presented  in  Plates  1  through  15.  In  each  plate, 

the  vertical  bar  connecting  H  ,  H,  ,  H~  ,  and  H  .  represents  the 

max  1/3  0.50  min  r 

physical  model  results,  and  the  circular  symbol  represents  the  numerical  model 
result.  In  each  case,  actual  wave  heights  have  been  normalized  by  the  height 
of  the  generated  wave,  as  listed  in  Table  1.  The  vertical  bar  plots  are  in¬ 
terpreted  in  Part  II.  Amplification  factors  from  the  numerical  model  were 
rounded  to  the  nearest  0.05  unit  of  the  incident  wave  height  before  multiply¬ 
ing  by  the  appropriate  shoaling  coefficient.  Most  of  the  comparisons  agree 
reasonably  well.  Surprisingly,  even  in  cases  such  as  Runs  2,  4,  8,  10,  and  14 
(corresponding  to  Photos  26,  38,  96,  102,  and  105,  respectively,  in  Bottin's 
(1984)  report)  where  wave  breaking  and  overtopping  occur,  the  comparisons  arc 
good  for  many  gages  though  they  are  poor  for  some.  These  poor  comparisons 
probably  result  from  the  effects  of  wave  breaking  and  overtopping  not  being 
strong  enough  in  these  cases  to  significantly  alter  major  features  of  the  har¬ 
bor's  response. 

42.  Most  of  the  numerical  model  wave  heights  are  between  H  and 

max 

Hq  from  the  physical  model.  However,  the  numerical  model  underpredicted 
wave  height  for  the  shortest-period  waves  (6.7-sec  waves  from  the  northeast  in 
Runs  5,  6,  11,  12,  and  15).  The  underpredictions  are  partly  caused  by  the 
fact  that  the  finite  element  grid  spacing  did  not  provide  adequate  resolution 
for  the  shortest-period  wave;  however,  for  the  7.5-sec  waves  from  the  North 
(Runs  3,  4,  9,  10  and  14)  the  comparisons  are  much  better  despite  the  fact 
that  the  resolution  was  only  slightly  better.  The  most  severe 


underpredictions  were  for  Runs  5  and  6,  which  were  made  using  the  base  Test  2 
harbor  configuration.  For  this  harbor  configuration,  waves  from  the  northeast 
can  easily  propagate  through  the  navigation  entrance  and  through  the  opening 
between  the  east  breakwater  and  the  city  dock.  Hence,  for  base  Test  2,  dif¬ 
ferences  in  the  direction  and/or  height  of  waves  incident  to  the  harbor  in  the 
physical  and  numerical  models  can  result  in  very  different  harbor  responses. 
For  Plan  42  and  Plan  58  differences  in  the  i.-.cident  waves  of  the  physical  and 
numerical  models  are  not  as  important  because  the  harbor  breakwaters  allow 
very  little  energy  into  the  harbor.  Wave  transformation  which  takes  place 
over  the  bathymetry  of  the  physical  model  between  the  wave  generator  and  the 
edge  of  the  finite  element  grid  is  not  simulated  by  the  numerical  model 
(paragraph  33) .  This  fact  may  partly  explain  the  poor  comparison  for  Runs  5 
and  6.  It  also  indicates  that  the  finite  element  grid  should  be  constructed 
to  include  as  much  of  the  important  bathymetry  as  possible. 

Irregular  Waves 


43.  For  the  irregular  wave  run,  the  incident  wave  spectrum  shown  in 
Fieure  20  was  used  as  input  to  the  numerical  model.  The  input  spectrum  was 
divided  into  5  frequency  bands  centered  at  0.116,  0.138,  0.161,  0.187,  and 
0.211  Hz.  At  each  gage,  the  spectral  density  at  these  five  frequencies  was 
predicted  by  the  numerical  model.  The  spectral  density  was  set  to  zero  at 
frequencies  below  0.09  Hz  because  the  incident  wave  spectral  density  was  zero 
at  lower  frequencies,  whence  the  HARBD  model  (which  is  based  on  the  linear 
theory)  must  yield  zero  spectral  density  below  0.09  Hz.  The  numerical  model 
study  did  not  consider  frequencies  greater  than  0.211  Hz  because  very  poor 
resolution  was  expected  using  the  finite  element  grids  shown  in  Figures  16 
through  18. 

44.  For  each  gage,  the  energy  spectra  yielded  by  the  physical  and  nu¬ 
merical  models  are  compared  in  Plates  16  through  29.  Notably,  the  scales  are 
not  the  same  for  each  spectrum.  The  comparisons  generally  agree  well  for  the 
range  of  frequencies  considered  by  the  numerical  model.  For  Gage  1  near  the 

harbor  entrance,  spectra  from  the  physical  and  numerical  models  are  very  simi- 

2 

lar  in  shape.  The  spectral  density  peaks  are  about  the  same  value  of  4.5  ft 

sec.  For  Gages  2  through  14  inside  the  harbor,  results  from  both  models  indi- 

2 

cate  that  spectral  densities  are  all  small  and  less  than  0.5  ft  sec  for 
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frequencies  greater  than  0.09  Hz.  This  occurrence  indicates  that  only  a  very 
small  portion  of  the  energy  of  the  incident  wave  train  reaches  inside  the 
harbor . 


45.  Differences  in  the  spectra  occur  at  frequencies  less  than  0.09  Hz. 
Below  this  frequency  the  numerical  model  spectra  must  predict  zero  spectral 
density  at  each  gage  because  spectral  density  of  the  input  spectrum  was  zero 
below  0.09  Hz  (paragraph  43).  However,  the  physical  model  spectra  show  sig¬ 
nificant  spectral  density  below  0.09  Hz,  particularly  in  the  very  low  fre¬ 
quency  end.  In  the  harbor,  the  wave  energy  at  lower  frequencies  may  be  asso¬ 
ciated  with  excitation  by  wave  groups,  nonlinear  wave  energy  transfer,  and 
breaking  and  overtopping. 
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PART  IV:  CONCLUSION 

46.  Most  of  the  wave  records  from  the  physical  model  tests  exhibited 

the  effects  of  nonlinearity  and  some  degree  of  wave  grouping.  In  several 

tests  wave  breaking  and  overtopping  occurred  which,  in  turn,  generated  low- 

amplitude,  high-frequency,  random  oscillations.  The  spread  of  the  wave  height 

distribution  at  many  gages  was  fairly  wide.  It  has,  therefore,  been  useful 

and  convenient  to  characterize  wave  height  at  each  gage  by  the  vertical  bar 

plot  of  H  ,  H  .  ,  H.  ,  and  H.  ,  rather  than  by  a  single  value 

max  min  1/3  0.50  7  6 

such  as  .  For  the  monochromatic  tests,  the  vertical  bar  plots  were  used 

to  compare  results  from  the  physical  and  numerical  models. 

47.  The  HARBD  numerical  model  used  in  this  study  is  based  on  linear 
wave  theory  which,  in  a  rigorously  theoretical  sense,  is  not  intended  to  model 
the  wave  conditions  of  the  physical  model  tests.  (Some  of  the  physical  model 
tests  include  nonlinear  wave  transformation  and  extensive  breaking  and  over¬ 
topping.)  Nevertheless,  results  of  the  two  models  are  reasonably  consistent 
in  that  they  yield  comparable  harbor  responses  for  most  cases  even  though  the 
magnitude  of  the  response  at  various  locations  in  the  harbor  was  very  differ¬ 
ent  in  a  few  instances.  These  differences  are  important  because  they  may 
materially  affect  the  ultimate  design  of  the  harbor.  It  must  be  emphasized 
that  the  comparisons  were  not  straightforward  because  the  physical  model  re¬ 
sults  are  characterized  by  a  distribution  of  wave  heights,  whereas  the  HARBD 
model  ’-esult  is  simply  the  height  of  a  monochromatic  wave. 

48.  For  the  monochromatic  wave  comparisons,  the  HARBD  numerical  model 

yielded  wave  heights  which,  in  most  cases,  were  between  H  and  H_  en 

max  0 . 50 

from  the  physical  model.  However,  for  the  shortest-period  waves — 6.7  sec  in 
Runs  5,  6,  11,  and  15 — the  wave  heights  from  the  numerical  model  were  less 
than  those  from  the  physical  model.  This  is  mainly  because  the  grid  spacing, 
which  was  approximately  one-fifth  of  the  incident  wave  length  for  these  waves, 
was  too  coarse  to  provide  adequate  resolution.  A  grid  spacing  less  than  one- 
seventh  of  the  incident  wave  length  is  recommended  to  obtain  the  best  results 
for  most  engineering  practice.  The  worst  predictions  were  for  Runs  5  and  6. 
The  large  discrepancies  between  wave  heights  from  the  physical  and  numerical 
models  for  these  runs  may  be  partly  due  to  wave  transformation  between  the 
wave  generator  pit  and  the  edge  of  the  finite  element  grid. 

49.  For  irregular  wave  comparisons  the  HARBD  model  yielded  spectra 
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which,  for  frequencies  greater  than  about  0.09  Hz,  agreed  reasonably  well  with 
those  from  the  physical  model.  However,  there  was  a  great  difference  for  the 
portion  of  the  spectra  below  0.09  Hz.  The  physical  model  showed  a  significant 
amount  of  spectral  density;  whereas  the  numerical  model  predicted  no  spectral 
density  since  there  was  zero  spectral  density  below  0.09  Hz  in  the  incident 
spectrum.  The  significant  amount  of  spectral  density  at  low  frequencies  may 
be  due  to  excitation  by  wave  groups,  nonlinear  wave  energy  transfer,  and  wave 
breaking  and  overtopping.  The  HARBD  numerical  model  is  incapable  of  simu¬ 
lating  these  phenomena. 

50.  Both  models  have  limitations  related  to  their  ability  to  accurately 
simulate  the  prototype.  For  the  physical  model  these  are  primarily  related 
the  scale  effects,  limited  basin  size,  performance  of  the  wave  generator,  and 
determination  of  reflective  properties  of  the  shore  and  harbor  structures. 

The  negative  effects  of  these  limitations  always  are  minimized  to  the  extent 
possible  as  was  the  'ase  during  the  physical  model  study  of  Barcelona  Harbor. 
The  main  limitation  on  the  HARBD  model  is  that  it  is  based  on  linear  wave 
theory  with  mild  bottom  slope  approximation.  Accuracy  of  the  HARBD  model 
decreases  as  nonlinearity  of  the  wave  increases.  HARBD  is  not  intended  to 
model  certain  phenomena  such  as  wave  breaking  and  overtopping.  Nevertheless, 
the  study  indicates  that  if  these  phenomena  are  not  dominant,  HARBD  is  able  to 
predict  wave  heights  inside  the  harbor  fairly  accurately.  For  engineering 
purposes,  HARBD  can  be  used  somewhat  beyond  the  limits  of  its  theoretical 
range  of  applicability.  For  some  problems  the  accuracy  of  the  solution  is 
adversely  affected  if  the  finite  element  grid  is  not  extended  sufficiently  far 
seaward  to  cover  the  entire  area  of  important  bathymetry.  The  error  is  not 
inherent  to  the  HARBD  model;  rather,  it  results  when  resources  available  to 
the  study  limit  the  size  of  the  finite  element  grid  which  can  be  constructed. 
This  problem  is  analogous  to  the  finite  basin  size  limitation  of  physical 
models . 

51.  HARBD  represents  a  significant  improvement  over  previous  numerical 
models  because  it  is  able  to  simulate  bottom  friction  and  boundary  absorption. 
HARBD  can  predict  wave  response  inside  a  harbor  more  accurately  than  previous 
numerical  models,  particularly  at  or  near  the  resonant  modes  of  the  harbor. 

52.  Since  no  prototype  data  were  available,  inferences  cannot  be  made 
concerning  the  ability  of  either  model  to  predict  wave  response  in  the  proto¬ 
type.  To  more  rigorously  assess  the  models,  prototype  data  are  indispensable. 


Furthermore,  conceiving  and  conducting  this  study  after  the  physical  model 
study  had  been  completed  presented  several  problems.  The  principal  problem 
was  that  physical  model  results  were  not  ideally  suited  for  the  purpose  of 
making  comparisons  with  results  from  the  HARBD  numerical  model.  For  example, 
some  waves  considered  in  the  physical  model  study  had  heights  and  periods 
which  indicated  they  might  not  be  well  described  by  the  linear  theory.  An 
intercomparison  study  similar  to  the  one  described  in  this  report,  but  with 
the  physical  and  numerical  models  studies  planned  and  conducted  simultane¬ 
ously,  would  be  useful.  Ideally,  prototype  data  would  also  be  available  for 
such  a  study. 

53.  Even  though  the  physical  model  results  were  not  ideally  suited  for 
making  comparisons  with  the  linear  HARBD  model,  some  differences  in  the  re¬ 
sults  of  the  two  models  stand  out.  Much  of  the  difference  in  the  monochro¬ 
matic  tests  can  be  attributed  to  the  long-period  beats  seen  in  the  physical 
model  wave  records.  Most  of  the  difference  in  the  irregular  wave  tests  is 
seen  at  the  extreme  low  frequency  end  of  the  spectrum,  below  the  low-frequency 
cutoff  of  the  generated  spectrum.  It  is  not  known  whether  these  are  nonlinear 
phenomena  which  also  occur  in  the  prototype.  If  they  are  spurious,  the  proce¬ 
dure  for  conducting  physical  model  studies  should  be  reassessed.  If  they  are 
not  spurious,  then  an  even  greater  emphasis  should  be  placed  on  developing  a 
nonlinear  numerical  model. 

54.  The  SPM  (1984)  wave  diffraction  method  requires  great  simplifica¬ 
tion  of  most  harbor  wave  problems.  In  the  case  of  Barcelona  Harbor  it  ex¬ 
hibited  serious  limitations,  particularly  in  neglecting  wave  reflections  from 
interior  harbor  boundaries.  Since  interior  reflections  were  the  primary 
source  of  problems  in  the  prototype  harbor,  the  SPM  method  is  considered 
unsuitable  for  this  study. 

55.  In  a  practical  engineering  project  study,  the  modeling  approach 
must  be  carefully  selected.  The  nature  of  the  important  hydrodynamic  pro¬ 
cesses  (e.g.  wave  reflection,  breaking,  overtopping,  and  resonance)  and  the 
scale  and  characteristics  of  the  harbor  area  are  critical  in  model  selection. 
Other  important  considerations  are  the  number  and  types  of  alternatives  to  be 
investigated,  the  importance  of  tangible  demonstrations  for  co-sponsors  and 
the  public,  and  the  available  time  and  funding.  The  primary  driving  force 
behind  the  choice  of  an  appropriate  modeling  approach  should  be  the  required 
accuracy  and  completeness  of  the  predictions.  Physical  models  can  yield 
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information  concerning  the  effects  of  wave  grouping,  nonlinearity,  wave  break¬ 
ing,  and  overtopping  which  are  not  available  from  numerical  models.  They  can 
also  provide  information  on  wave-induced  currents  and  shoaling  patterns.  For 
some  studies,  numerical  models  have  advantages  in  terms  of  time,  cost  and 
flexibility.  SPM  (1984)  methodologies  are  useful  for  rough  estimates  and  for 
very  quick  studies  of  simple  harbors  but  have  serious  limitations  for  most 
actual  project  studies.  This  report  focuses  on  the  accuracy  of  the  various 
modeling  approaches.  It  does  not  provide  detailed  information  on  the  time  and 
cost  requirements  of  physical  and  numerical  models. 
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APPENDIX  A:  ESTIMATION  OF  WAVE  CONDITIONS  IN 
BARCELONA  HARBOR  DSTNG  DIFFRACTION  DIAGRAMS 


1.  Diffraction  diagrams  presented  in  the  Shore  Protection  Manual*  (SPM) 
(1984)  were  used  to  estimate  wave  conditions  in  Barcelona  Harbor  for  Runs  1, 

3,  and  5.  (See  Table  1  (main  text)  for  characteristics  of  the  generated  waves 
of  these  runs.)  The  harbor  configuration  for  all  these  runs  was  the  base 
Test  2  harbor  configuration.  Navigation  entrances  of  the  Plan  42  and  Plan  58 
harbor  configurations  deviated  so  much  from  the  simple  gap  entrances  for  which 
the  diffraction  diagrams  were  developed  that  it  was  impossible  to  even  attempt 
to  use  diffraction  diagrams  to  estimate  wave  conditions  for  these  harbor  con¬ 
figurations.  It  was  realized  that  the  methods  in  the  SPM  (1984)  were  not 
applicable  even  for  the  Base  Test  2  harbor  configuration  because  most  of  the 
assumptions  upon  which  the  methods  are  based  were  not  satisfied  in  the  proto¬ 
type.  Nevertheless,  it  was  deemed  worthwhile  to  try  to  make  crude  estimates 
and  to  compare  the  results  with  results  from  the  numerical  and  physical 
models . 

2.  The  method  used  to  determine  wave  heights  inside  the  harbor  using 
diffraction  diagrams  follows.  First,  a  cartesian  coordinate  system,  at  the 
prototype  scale,  was  established  for  the  harbor,  the  origin  of  which  was  lo¬ 
cated  at  the  center  of  the  breakwater  gap.  The  y-axis  was  established  co- 
linear  with  the  bisector  of  the  angle  formed  by  the  two  breakwaters  with  y 
increasing  toward  shore.  The  orientation  of  the  x-axis  is  immaterial  since 
the  diffraction  diagrams  are  symmetric  about  the  y-axis.  Next,  the  coordi¬ 
nates  of  Gages  2  through  14  were  determined  (Gage  1  was  not  considered  since 
it  is  outside  the  harbor),  in  prototype  units,  for  the  stated  coordinate  sys¬ 
tem.  Then  for  Runs  1,  3,  and  5,  the  following  steps  were  used  to  determine 
the  wave  height  at  Gages  2  through  14: 

a.  The  linear  theory  dispersion  relation  was  used  to  calculate  the 
wavelength  using  the  period  of  the  generated  wave  and  the  depth 
at  the  breakwater  gap. 

b.  The  coordinates  of  the  gages  were  nondimensionalized  by  dividing 
by  the  wavelength  from  a. 


Referencp«  rit-ed  in  the  Appendix  can  be  found  in  the  References  at  the  end 
of  the  main  text. 


c.  The  ratio  of  gap  width  to  wavelength  B/L  was  determined,  and 
this  value,  along  with  the  direction  of  the  generated  wave,  was 
used  to  select  the  appropriate  diffraction  diagram  in  the  SPM 
(1984). 

d.  The  selected  diffraction  diagram  along  with  the  dimensionless 
gage  coordinates  were  used  to  determine  the  diffraction  coeffi¬ 
cient  for  each  gage. 

e.  The  wave  height  at  each  gage  was  calculated  as  the  product  of 
the  diffraction  coefficient  for  that  gage  times  the  incident 
wave  height  times  a  shoaling  coefficient  (see  paragraph  33,  main 
text) ;  the  latter  accounted  for  shoaling  between  the  wave  gen¬ 
erator  and  breakwater  gap.  Shoaling  coefficients  ranged  from 
1.03  to  1.14. 

3.  Wave  heights  thus  determined  are  listed,  along  with  values  found 
from  the  physical  and  numerical  models,  in  Tables  A1 ,  A2,  and  A3  for  Runs  1, 

3,  and  5,  respectively.  The  SPM  method  includes  no  reflection  and  hence  sig¬ 
nificantly  underestimates  wave  heights  in  some  areas,  particularly  at  the  city 
dock  (Gages  11  and  12). 


Numerical  Model,  and  Physical  Model* 
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Values  given  are  predicted  wave  height  for  the  physical  model) 

normalized  by  the  generated  wave  height  from  Table  I  (main  text) . 


Table  A2 

Comparison  of  Wave  Heights  for  Run  2  Using  Diffraction  Diagram, 
Numerical  Model,  and  Physical  Model 
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Table  A3 


Comparison  of  Wave  Heights  for  Run  3  Using  Diffraction  Diagram, 


Numerical  Model,  and  Physical  Model 


